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EDGE STATISTICS FOR A CLASS OF REPULSIVE PARTICLE 

SYSTEMS 

THOMAS KRIECHERBAUER AND MARTIN VENKER 


Abstract. We study a class of interacting particle systems on M which was recently inves¬ 
tigated by F. Gotze and the second author |GV14| . These ensembles generalize eigenvalue 
ensembles of Hermitian random matrices by allowing different interactions between parti¬ 
cles. Although these ensembles are not known to be determinantal one can use the stochastic 
linearization method of | GV14| to represent them as averages of determinantal ones. Our 
results describe the transition between universal behavior in the regime of the Tracy-Widom 
law and non-universal behavior for large deviations of the rightmost particle. Moreover, 
a detailed analysis of the transition that occurs in the regime of moderate deviations is 
provided. We also compare our results with the corresponding ones obtained recently for 
determinantal ensembles |Schl5llEKS15| . In particular, we discuss how the averaging effects 
the leading order behavior in the regime of large deviations. In the analyis of the averaging 
procedure we use detailed asymptotic information on the behavior of Christoffel-Darboux 
kernels that is uniform for perturbative families of weights. Such results have been provided 
by K. Schubert, K. Schuler and the authors in | KSSV14] . 


1. Introduction and Main Results 

In [BdMPS95] . a class of particle ensembles on the real line with many-body interac¬ 
tions was introduced, generalizing invariant random matrix models, and results on the global 
asymptotics were obtained. Gotze and the second author |GV14| presented a stochastic lin¬ 
earization procedure which allowed to study local correlations of particles in the special case 
of two-body interactions. They showed that in the limit of infinitely many particles the bulk 
correlations are described by the celebrated sine kernel law. The purpose of this paper is to 
extend the investigations of |GV14| to the edge. 

To be more precise, we consider particle ensembles on M with density proportional to 

( 1 ) 

i<j 

where Q : R —>• R is a continuous function of sufficient growth at infinity compared to the 
continuous function (p : R —)• [0,oo). Apart from some technical conditions we will assume 
that 

= 0, pit) > 0 for f 7 ^ 0 and lim = c > 0, 

^ ^ |t|2 
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or, in other words, 0 is the only zero of and it is of order 2. The connection to random 
matrix theory is as follows: If ^p{t) = then we have the unitary invariant ensembles Pn,q 
given by 


Zn.Q y 

Here and in the following, we slightly abuse notation by using the same symbols for a measure 
and its density. Pn,q is the joint distribution of the eigenvalues for the random (N x N) 
Hermitian matrix with density proportional to 

g-VTr(Q(A)) 

w.r.t. the Lebesgue measure on the space of {N x N) Hermitian matrices. Here Tr denotes 
the trace and Q{X) is defined by spectral calculus. This matrix distribution is invariant 
under unitary conjugations. See e.g. [Dei99[[ABDFllllPSlllE.GZ10] for general references on 
random matrix theory. The unitary invariant ensembles ([3]) have been well-understood in 
terms of asymptotic behavior of the eigenvalues. A common way to describe their properties 
is to use correlation functions. For a probability measure P]\f{x)dx on and a natural 
number k < N define 

p%{xi,... ,Xk) := / PN{x)dxk+i ..-dxN- 

The k-th correlation function is a probability density on M^, the corresponding measure is 
called fc-th correlation measure. Note that in the Random Matrix literature the fc-th correla¬ 
tion function often differs from the definition above by the combinatorial factor A^!/(A^ — k)\. 
The key for the detailed understanding of the local eigenvalue statistics of unitary invariant 
ensembles Pn,q is that they are determinantal ensembles with a kernel that can be expressed 
in terms of orthogonal polynomials. More precisely, denote by the Christoffel-Darboux 

kernel of degree N associated to the measure e~^^^^'^dx, i.e. 

N-l 

where is the unique polynomial of degree j with positive leading coefficient such that 

forms an orthonormal sequence in L‘^{e~^^^^^dx). Then, for all 1 < A: < one 
can express the k-th. correlation function by 

A^! k 

_j^y PN,Q{xi,...,Xk) = det {{KN,Q{Xi,Xj))i<ij<k) . (4) 

To state our results, it is convenient to rewrite ©• Let h : M —M be continuous, even and 
bounded below. Let Q : M — > M be continuous, even and of sufficient growth at infinity. 
Consider the probability density on given by 

1 ^ 

Pn,q{x) := -p—-a;j|%xp{-A^^Q(xj) - ^/i(xi - Xj)}, (5) 

^iQ i<j j=l i<j 
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where q is the normalizing constant. Choosing ip{t) := exp{—we see that q 
is in the form ([T]). Conditions ([2]) are satisfied as h has no singularities. 

G. Borot has pointed out to the second author that q is the eigenvalue distribution of 
the unitary invariant ensemble of Hermitian matrices with density proportional to 

^-N'Tr{Q{X))-Tr{h{X®I-I®X)) 

where I denotes the identity matrix (cf. [BE015] i. Nevertheless, we will prefer to view Pj^ q 
as particle ensemble and therefore speak of particles instead of eigenvalues. Furthermore, the 
ensemble does not seem to be determinantal, in contrast to Pn,q- 

Our method of proof requires the external field Q to be sufficiently convex. To quantify 
this, we define for Q being twice differentiable and convex ag := inf^gR Q"{t). In the following 
we will denote by the /c-th correlation function or measure of P^ q . 

The first information we need is the global behavior of the particles, i.e. their asymptotic 
location and density. To this end, we state the following theorem which was proved in |GV14] . 

Theorem 1 ( |GV14t Theorem 1.1]). Let h be a real analytic and even Schwartz function. 
Then there exists a constant > 0 such that for all real analytic, strictly convex and even 
Q with aq > , the following holds: 

There exists a compactly supported probability measure pq having a non-zero and continuous 
density on the interior of its support and for k = 1,2,..., the k-th correlation measure of 
P^ q converges weakly to the k-fold product of pq, that is for any bounded and continuous 
function g : —> M, 

J gih, ■ ■ ■ ,tk) PN%iP^ ■ ■ ■ Pk) d'^t = j g{ti,... ,tk) dpqih)... dpqitk). ( 6 ) 

Remark 2. 

a) If h is positive-definite (i.e. with nonnegative Fourier transform), can be chosen as 

= supjgR —/i"(t) (cf. Remark [14] and |GV14[ Remark after Theorem 1.1]). 

b) It was shown in |GV14l Sec. 3] that the measure pq is the equilibrium measure with 
respect to the external field V{x) := Q(x) + f h{x — t)dpq{t). This fact leads to a 
useful representation of pq (see e.g. AppendixjAj). We note for later reference that the 
evenness of Q and h carries over to pq (cf. |GV14[ Theorem 1.1]) and hence to V. 

c) In |BdMPS95| . ensembles with many-body interactions are considered and global 
asymptotics but not local correlations are derived (see also the discussion on page 
326 in |TotOO| f . In the case of pair interactions, the classes of admissible interactions 
in |BdMPS9^ and in this paper are different. In |BdMPS9^ . a convexity condition is 
posed, depending solely on the additional interaction potentials where our conditions 
depend on both Q and h. The characterization of the limiting measure is different, 
too. 

Asymptotic expansions of the partition functions can be found in |BGK13] . Further 
global asymptotics for repulsive particle systems have been studied in [GGZ141IGP14] . 

On the local scale, there are basically two different regimes of interest. The local regime 
in the bulk of the particles has been investigated in |GV14j . yielding asymptotics in terms 
of the sine kernel. These results on the bulk have been extended in two ways. Firstly, 
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in the recent [SV15j the sine kernel asymptotics have been obtained under an unfolding of 
the particles, which is a non-linear rescaling that transforms the equilibrium measure to a 
uniform distribution. This in particular allows for an efficient treatment of empirical spacings 
between particles. Secondly, bulk universality has been shown in [Venl3| for the /3-variant of 
(jS]) which is obtained by replacing the exponent 2 in ([5]) by some arbitrary /3 > 0. 

Returning to /3 = 2, the asymptotics at the edge are governed in terms of the Airy 
kernel. The Airy function Ai : M —M is uniquely determined by the requirements 
that it is a solution to the differential equation f"{t) = tf{t) and has the asymptotics 
Ai(t) (47r\/t) 3^^^^ as t —>• oo. The Airy kernel ATai is defined on and can be 

represented by 


KAiit,s) = 


POO 

/ Ai(t 

Jo 


+ r)Ai(s -|- r)dr = 


Ai(t)Ai'(s) — Ai'(t)Ai(s) 
t — s 


( 7 ) 


where the latter representation only holds for t ^ s. The first main result of the present work 
is the edge universality of the correlation functions of q- Here the universal aspect is that 
in the limit N ^ oo the correlations do not depend on Q and h after an appropriate linear 
transformation t b + t/{c*N‘^P). The quantities b and c* depend on Q and h through the 
function V defined in Remark [2] b) as follows. The number b is the supremum of the support 
of the equilibrium measure fig and c* is related to the prefactor of the square root vanishing 
of the density of Hq at 6, 


9 - 1/3 

C* := ——Gy(l)2/3, 

b 

where the function Gy is defined in (I75p . 


( 8 ) 


Theorem 3. Let h and Q satisfy the assumptions of Theorem [7] and assume furthermore 
that the Fourier transform h of h decays exponentially fast. Let q gM be given and let b and 
c* be as above. Then we have for N ^ oo 

a) 

Pn% (b + ...,b + = det [KAiiti, + o(l) (9) 

with the o(l) term being uniform in t G [q,oo)^. 

b) If h is negative-definite, i.e. h<0, then for any 0 < a < 1/3 the o(l) term in Q can 
be replaced by 0{N~^), again uniformly in t G [g, oo)^. 


Remark 4. 

a) Even in the situation of part b) where we have better bounds on the rate of convergence, 

the Airy kernel determinant is the leading term in ([9]) only for t G [q, o{{log . 

b) The condition of exponential decay of the Fourier transform of h allows for writing h 
as difference of positive-definite and real-analytic functions. This technical condition 
is needed and explained in Section [3l In particular, if h is positive-definite or negative- 
definite, the assumption on the decay of the Fourier transform is automatically satisfied 
[LS52]. 
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c) If in the situation of part a) of the previous theorem, h is positive-definite, then 
can be chosen as = sup^gR —(cf. Remark [2]). This remark extends to all of 
our results. The roles of the different conditions on the dehniteness of h are explained 
in Remark 1141 below. 

d) Edge universality has been established for a wide variety of random matrix ensembles. 
For /3-variants of ([3]) with one-cut support we refer the reader to the recent works 
|KRV15| . [BEY14] and |BFG15| . It should be noted that the limit laws depend on the 
value of /? [RRVll] . For results on invariant matrix ensembles and on Wigner matrices 
we refer to the books and recent surveys [DGOhllPSlllITVldllErdM] and references 
therein. 


We now turn to the largest particle Xmax and provide an extensive description of its dis¬ 
tributional properties. In analogy to the classical central limit theorem for independent 
variables, we prove a limit theorem on the fluctuations of Xmax around the rightmost end¬ 
point b (Theorem [5|). It also implies the weak law of large numbers with respect to a natural 
product measure. The strong law of large numbers is also shown (Corollary [TT|) . Special 
emphasis is given to the study of the upper tail of the distribution function of Xmax- Strong 
results on moderate and large deviations are presented in Theorem [6] and Corollary [HI We 
already point out, that it is only in the regime of large deviations that the differences between 
the repulsive particle systems q and the determinantal ensembles become visible in the 
leading order behavior. 

We start with the fluctuations of Xmax around b. The limiting law is not Gaussian but 
is given by the Tracy-Widom distribution. Its distribution function F 2 can be defined as a 
Fredholm determinant implying in particular the representation 


F2{s) 





It can also be expressed as 


F2{s) 



{t - s)q{tfdt}, 


( 10 ) 


where q is the solution of the differential equation q"{t) = tq{t) + 2g(t)^ which has the 
asymptotics q{t) ~ Ai(t) as t —>■ 00 |TW94| . For later use we note that the large-s-asymptotics 
of F 2 {s) are known (see e.g. [BBDOSl (I), (25)]) to be 

Theorem 5. Let Xmax denote the largest component of a vector x. Under the conditions and 
with the notation of Theoreml^ we have for iV —>■ 00 

a) 


Pn,q ((a^max - b) c*N^/^ <s)= F 2 {s) + 0 ( 1 ) 


( 12 ) 


uniformly for s E M, where F 2 is the distribution function of the ((3 = 2) Tracy-Widom 
distribution. 
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b) If h is negative-definite, then the o(l) term in (fT^ can be replaced by 0{N ^), where 
0 < (T < 1/3. The O term is uniform for s E [q, oo) with arbitrary gr E M. 

Note that Theorem [5] implies the weak law of large numbers w.r.t. the product measure 
IP •= 0A=i Q ™ appropriate measurable space), i.e. Xmax converges in probability to 
b. The strong law, i.e. almost sure convergence, is the content of Corollary 1111 below. 

Theorem [5] can be viewed as an analogue of the central limit theorem for the distribution 
of the largest particle. This point of view leads naturally to the question of what can be 
said about its large and moderate deviations? There one is interested, roughly speaking, 
in describing the decay of the probabilities P^q ((xmax — b) c*N‘^/‘^ < —s) (lower tail) and 
Q ((^max “ b) c*N'^P > s) (upper tail) when s is allowed to grow to infinity at some N- 
dependent rate. In this paper we deal exclusively with the upper tail. Note that even for 
determinantal ensembles the lower tail has not been studied at the level of detail required for 
our analysis. 

The large (and moderate) deviations theory generally deals with probabilities decaying at 
an exponential rate. In a general setting with probability measure P and random variables 
X ]\[, the deviation s is often parametrized as s = rfor some r, a > 0 and one proves 
statements like 


lim = -/(.). (13) 

N^OO ^ ^ ^ ^ 

Then N'^ is called the speed and I{r) the rate function. In our situation P corresponds 
to P ;= (S)'n=i q corresponds to (xmax — b) c*N‘^P with x E The regime of 

moderate deviations extends over 0 < a < 2/3 and the regime of large deviations is associated 
with a = 2/3. Note that Theorem [5] does not provide any information on either the speed or 
on the rate function for any a > 0 due to the fast decay of 1 — F 2 described in (fTTI) . 

The next theorem and its two subsequent corollaries remedy this situation. E.g. in the 
case of moderate deviations Corollary [TO] b) shows for all 0 < a < 2/3 that we have speed 
and rate function I(r) = For large deviations (a = 2/3) we learn from Corollary 

m a) that the speed is N and the rate function /(r) = r]v{l + r/{bc*)) with tjv defined in 
m- Observe that in the moderate regime the rate function is universal, whereas the rate 
function depends on V and is hence non-universal for large deviations. 

The situation becomes more subtle, if one studies the asymptotics of the tail probabilities 
and not just their logarithms. These provide finer information about the tails and we show 
how the E-dependence of the leading order behavior becomes stronger within the region of 
moderate deviations (Remark [9|). 

In recent years superlarge deviations have been studied as well, which would correspond 
to a > 2/3. We will not treat them here. See Remark [7] for a description of the difficulties 
that appear in that regime. 

We now state Theorem El that summarizes our results on moderate and large deviations 
for the upper tail of the largest particle. It is stated for Pj!)Q(xmax > t) and hence s in 
the discussion above corresponds to t = b -h s/{c* In this scaling the regime of large 
deviations is described (more naturally as for s) by those values of t for which f — 6 is of order 
1 and the regime of moderate deviations reads N~‘^P « t — b « 1. Theorem [6] provides 
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upper and lower bounds on Q{xina,x > t) for t belonging to an arbitrary but fixed bounded 
subset of (6 + oo). 

In the regime of moderate deviations the quotient of the npper and lower bound is seen to 
converge to 1 as —)• oo and thus the leading order behavior of -P^ ^(xmax > t) is determined. 
It is given by 

y2^-Nr]v{t/b) 

•= ATrN{P - b^)r,'y{t/b)' 

The ^-dependence is contained in the function rjv ■ [l,oo) —)■ M, 


riv{t) ■■= / Gvis)ds, 


(15) 


where the function Gy is defined in (j75[) below for the external field V introduced in Remark 
[2] above. In the regime of large deviations, onr resnlt is weaker. There we show that the 
quotient of pN,v{t) and q{x max > t) and its inverse are both bonnded. 

Theorem 6. Assume that the conditions of Theoreml^ are satisfied and fix T > b. Then the 
following holds: 

a) There exists a constant B > 1 such that for all N and all t ^ {b + N~'^/^,T) we have 

< PN,Q{^raax > t) < BF^yit)- (16) 

b) For any sequences {pn)n, {<1n)n of reals with pjy < qn, N‘^^^{pn — b) ^ oo, and 
qn — b ^ 0 for N ^ oo we have 

^N,Q (^^max > t) = -pAf,u(l)[l + o(l)] (17) 

uniformly in t G {pn-, Qn) as N ^ oo. 

c) If h is negative-definite, then o(l) in statement b) can be replaced by the more precise 


O 


1 


+ 0[y/t-b 


N\t-b\^/‘^ ^ 

with the O terms being uniform in N and in t G {bN~‘^^^,T). 


(18) 


Note that for negative-definite functions h statement c) implies the upper bound in (I16p 
but not the lower bound. 


Remark 7. Theorem [6] should be compared with the result |Schl5l (1.14)] (cf. |EKS15j l for 
determinantal ensembles. Even our best error bounds in the case of negative-definite functions 
h are worse than the bonnds there, if \t — b\ » . This is due to the 0{y/t — b) term 

in ([18]) . We believe that this is not an artefact of onr proof. In fact, relations (1621) and (l63]) 
of the proof show that we can express the distribntion of the largest particle by averaging 
determinantal ensembles with V replaced by V — f /N over /. Tracing back the origin of the 
C>(Vl — b) term in (I18jl we see that it does not stem from a lack of control on error terms but 
on the dependence of the upper endpoint oa /, that can actually be computed to first 

order. In the regime of large deviations this effect is most prominent and this is the reason 
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why we loose the convergence result for the distribution of the largest particle. However, we 
conjecture that also in the regime of large deviations we have convergence 

lim Pko (s^max ^ — P (J'')Pn,V (i) 

N—>-oo 

for some strictly monotone increasing function B with B{t) > 1. Observe that our previous 
results all show that q has the same limiting behavior as the determinantal ensemble Pn,v- 
However, for Pn,v one has B{t) = 1, i.e. we conjecture the breakdown of comparability 
between Pj^^Q and Pn,v in this regime! To show this one would need to apply Laplace’s 
method to the infinite dimensional integral that represents the averaging over /. This is a 
technically demanding enterprise but not without hope due to the Gaussian character of the 
probability measure on /. We plan to come back to this problem in future work. 

The situation becomes even more challenging in the superlarge regime, where the additional 
difficulty arises that one may not truncate Pj!f q to some finite interval [—L, L]^ that is used 
in the proof of Theorem [6] to obtain good bounds on the Gaussian field / used for the 
linearization method (see e.g. discussion in the paragraph before (1461) L In contrast, for 
determinantal ensembles the leading order behavior of the upper tail Pn,v (a^max > t) was 
shown to be given by pN,v{t) in both the large and the superlarge regime for a suitable class 
of convex functions V |Schl5( Theorem 1.1] (cf. [EKSlhj L 


Next we discuss the connection with the Tracy-Widom law F 2 beyond the central limit 
regime already covered in Theorem [5l It is convenient to use the scaling of Xmax as in the 
statement of Theorem (Sj i.e. t = b + s/(c*N‘^^^). The positivity and real analyticity of 
Gv (see ([7^ and thereafter), definitions ([TH]) . (fT4j) and the definition of c* in ([8]) lead by a 
straightforward calculation to the following representation of Pn,v in the regime of moderate 
deviations, i.e = o(l): 


pN,v{t{s)) 

Nrjv{t{s)/b) 


Q-Nriv{t(s)/b) 

167rs3/2 + ^ 



CX> 


+ X] 

i=i 


s^'+f 



with 


= -s3/2 + 0 
3 



(19) 


for some sequence {djy)j>i of real numbers depending on V. 
thank F. Gotze for pointing out that the power series 



00 


+ X] 

i=i 




The authors would like to 


is an analogue to the Gramer series |Cra38] in the deviations theory for sums of independent 
random variables. It follows that 

e-N^vm/b) =e-|*'^'[l + o(l)] 


as —>■ 00 for 0 < s = o (iV^/^®), or equivalently for 0 < t{s) — b = o [N Gomparing 

with the asymptotics of the Tracy-Widom distribution (Hip and using in addition that 



for s > 0 
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we obtain from statements b) and c) of Theorem [6l 

Corollary 8. Under the condition and with the notation of Theorem\^ the following relations 
hold. 

a) For any sequences {pn)n, iQN)N of reals with < q^, p^ oo, and N~^/^^q]\[ 0 

for N ^ oo we have: 


Pi 


N,Q 


((W - 6)c*iV2/3 > s) = (1 - F2 (s))( 1 + o(l)), 


uniformly in s £ {pN, qw) as N ^ oo. 
b) If h is negative-definite, then 


Pn,q ((w - 6)c*iV2/3 > s) = (1 - F2{s)) 1 + O 
with the O terms being uniform in N and in s £ 


+ 0 


g5/2 

Ar2/3 


Remark 9. With Corollary [5] we have identified the region for which the leading order of 
the upper tail of the distribution of the largest particle is universal, i.e. all of its Q and 
h dependence is encoded in the two numbers b and c* that define the linear rescaling and 
play the same role as mean and variance in the Central Limit Theorem. Relations (I19p show 
nicely how the dependence of the leading order pNy{t) on Q and h grows gradually within 
the regime of moderate deviations. Set 7 ^ := 2l{2k + 5). For 0 < t — b = i.e. 

s{t) = o (A^“'''fe+2/3^ ^ function pNyif) depends to leading order on Q and h only through 
the k + 2 numbers b, c*, diy,... ,dky. This is a generic phenomenon already observed 
in [Schl5[ Remark 4.12] (cf. [EKS15] ) for determinantal ensembles. 


In the logarithmic form, it is straightforward to derive the following results from Theorem 
[6] and relations (fT^ , (fTTl) , (fT^ and (1751) . 


Corollary 10. Under the condition and with the notation of Theoreml^ the following state¬ 
ments hold. 


a) 




N 




where the O term is uniform in N and in t £ (6 + N~'^P,T). 
b) For any sequences {pn)n, iQN)N of reals with pN < qN, Pn —t 00 , and N~‘^Pq]s[ 0 

for N ^ 00 we have 


logP^ Q ((Xmax - b)c*N‘^P > s) 


_4 _ log(167rs^/^) / J_\ / s \ 

3 s3/2 ^ Vs3/2y 


with error bounds that are uniform in N and s £ {pn,Qn)- If h is negative-definite, 
then the o (s“3/2^ term can be replaced by O that is also uniform in s £ {pN, §jv). 


Finally we present the strong law of large numbers for the largest particle that will be 
proved at the end of Section [H 
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Corollary 11. Assume the conditions of Theorem\^ and set P := Q- Then 

P( lim Xmax = 6) = 1. 

N^oo 

Remark 12. Corollary [TO] a) implies a full large deviations principle on M for Xmax with 
speed N and good rate function t i—>• r]v{t/b), where we set r]v{t) = oo for t < 1. Such a 
large deviations principle has already been proved in a more general setting for /3-variants 
of ([3]) (see [BGKl.lj and references therein). It can be derived via a contraction principle 
from a large deviations principle for the empirical distribution of the particles. We also refer 
the reader to [lrTo] and for results on combinatorial models that are related to invariant 
matrix ensembles to |Joh00[lBDM~*~0l[|LM01llLMR02] . As mentioned above, stronger results 
in non-logarithmic form as in Theorem [6| and Corollary [8| have been shown for determinantal 
ensembles in |Schl5[[EKS15| . 

The paper is organized as follows. Section [2] provides a sketch of the central ideas introduced 
in [GV14] that guide the proofs of all our main results. Then, in Section [S] we derive our 
universality result on the correlation functions. The last section deals with our various results 
on the distribution of the largest particle. In Appendixj^ we collect all results from [KSSVI4] 
on determinantal ensembles that are needed in this paper and bring them in a form that is 
suitable for our purposes. 


2. Outline of the Method 


We first mention the basic steps from |GVI4j in the analysis of P^ q- The additional 
interaction term h{xi — Xj) is in general of order so it may influence the limiting 

measure. The idea is to split it into a term contributing to the formation of the limiting 
measure and a perturbation term of lower order. To this end, let us introduce for a probability 
measure /r on M the notation := f h(t — s)dfi{s) and •= / / “ s)dfj,{t)dfi{s). We 

can then write 

h{xi - Xj) = \Y 
i<j id 

=N Y + 2 ^ ~ ~ hfj,{xi) - h^{xj) + + Cn, ( 20 ) 

3=1 id 

where Cn := —{N/2)h{0) — {N‘^/2)h^^. The term in brackets in (I^UD is the Hoeffding 
decomposition of the statistic h{xi — Xj) w.r.t. the measure /r. The term N V(^3) 

will be added to the external field Q forming a new potential ■.= Q + h^. Setting 

U^{x) [Hxi - Xj) - h^{xi) - hfj,{xj) + , ( 21 ) 

id 


we arrive at the representation 


PN,v,ix)e^-^P 


(22) 
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with •= ■ Formula (l2^ establishes a relation to a determinantal ensemble, 

indicating that asymptotics of q may be deduced from asymptotics of Pn,v^ ■ Our aim is 
to find a measure /r such that the ratio and its reciprocal are bounded in N. 

It is straightforward to see that 




^N,v^e' 




(23) 


where 1 Ea,V'^ denotes expectation w.r.t. Pn,v^- In view of the desired boundedness of (1251) . 
should be centered at least asymptotically, which motivates the condition that jjL should be the 
equilibrium measure to the field V)j. This implicit problem was solved in |GV14[ Lemma 3.1] 
by a fixed point argument, yielding existence but not uniqueness of such a /r. The uniqueness 
followed later by proving that any measure /r which is the equilibrium measure to is the 
limiting measure for P^ q- From now on let /i denote the unique measure with this property, 
write V := and U Note that this definition of V is consistent with that in Remark 

Elb). 

At this stage, it can be proved that the ratio Z^y/Zjyyu is bounded in N and bounded 
away from 0 provided that aq is large enough. More precisely, given A > 0, there is an 
a(A) < oo such that there are constants 0 < C'i(A) < C 2 (A) < oo such that for oq > a(A) 


Ci{X)<ENye^ < C2(A) 


(24) 


for all N (see |GV14[ Proposition 4.7] and [GV141 Remark 4.8]). 

One main tool to derive bound (I24p is the following representation of U in terms of linear 
statistics (cf. |GV14[ Lemma 4.6]) using Fourier techniques, 


U{x) 



2 




h{t)dt, 


(25) 


with h{t) = (27r) f e ^*^h(s)ds. The other ingredient to the proof of ([Mil is the following 
concentration of measure inequality for linear statistics which will be needed later on, too. 


Proposition 13. Let Q be a real analytic external field with Q" > c > 0. Then there exists 
a positive constant C such that for any Lipschitz function g whose third derivative is bounded 
on an open interval / C M containing the support of the equilibrium measure pg, we have for 
any e > 0 and N 


EN,Qexp{£{'^g{xj) - N j g{t)dpQ{t))] < exp{ ^ +eC(||g||oo + ||g^^^||oo)} 
i=i 


where \g\i^ denotes the Lipschitz constant of g and || • ||oo denotes the sup norm on L. 

The proposition can be found in [GV14t Corollary 4.4] and follows from the fact that for 
strongly convex Q, Pn,q fulfills a log-Sobolev inequality which implies concentration of the 
Lipschitz function around its expectation (see e.g. [AGZlOi Proposition 4.4.26]). 

To obtain Proposition 1131 we combine this with the following estimate on the distance of the 
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expectation to its large N limit from [Shell! Theorem 1] (see also [KS10[ Theorem 1]), 

N 

^N,Q'^g{xj) - N 
i=i 

where C is the constant that then appears in Proposition 1131 

The main idea to tackle the local universality is linearizing the bivariate statistics U to 
transform it into linear statistics. To this end, assume that —h is positive-definite (i.e. h < 0; 
the general case is reduced to this case). It is well-known (see e.g. [AT07| 1 that due to the 
positive-definiteness of —h, there exists a centered, stationary Gaussian process (/(t))teR 
which has —h as its covariance function. That is, the finite-dimensional distributions of / 
are Gaussian, E/(t) = 0 and E/(t)/(s) = —h{t — s), where E denotes expectation w.r.t. 
the probability space underlying the Gaussian process. The key observation now is that for 
fixed X E the random variable f{xj) ~ ^ J fdlJ- is again a centered Gaussian with 

variance 2U{x), yielding 

N „ N 

exp{7/(x)} = Eexp{^/(xj) -N fd^} = Eexp{^/(xj)} (26) 

i=i i=i 

with f := f - j fdfj. (27) 


gdg < Cdiffiloo + 


Gombining (j26l) with (|22l) and (1231) . the ensemble q can now be represented as an average 
over determinantal ensembles P/v,y-//iv a-s follows 


Pknix) — 


^ [^N,V-f/NPN,V-f/Nix)] 


EZ 


Ny-f/N 


(28) 


Relation (1281) immediately extends to correlation functions. At this point the key observation 
for our universality proofs is the fact that for almost all / the correlation functions v-f/N 
and Y coincide in the region of interest to leading order if N is sufficiently large. We will 
show in the appendix how to extract this information from [KSSV14] (cf. Proposition I16p . 

We end the brief overview of the method of proof by noting that in order to use the 
rich theory on Gaussian processes, in particular the sub-Gaussianity of maxima of Gaussian 
processes on a compact, we first need to truncate the ensemble to some compact interval 
[-L,L]^. 

Finally, we comment on the role that positive/negative-definiteness of h plays in the ap¬ 
proach of |GV14] and how it extends to the present paper. 


Remark 14. The reader might have wondered about the different results depending on h being 
positive-definite or negative-definite (cf. Remark |4|). If h is positive-definite, then 7/ in (1211) 
is negative for all configurations x, which can be easily seen from the Fourier representation 
([25]) . This immediately gives the upper bound of (IMl) . without any condition on ag. The 
lower bound and the further analysis in the case of positive-definite h only need V to be 
strongly convex, e.g. to apply Proposition [131 With the crude bound ag > sup^gj^ —/i"(t), 
given in Remark [2l this is guaranteed. 

If h is not negative-definite, then the linearization method as outlined in (|26[) needs to be 
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modified. The function h can be decomposed into positive-definite functions h^, such that 
h = h'^ — h~ . If 7 ^ 0, then —h is not a covariance function, but all functions zh'^ + h~ with 
z > 0 are and can be used for the linearization of a modified version Uz olU. We will then 
employ Vitali’s Theorem from complex analysis (a consequence of Montel’s Theorem and the 
Identity Principle) to transfer the convergence results for 2;>0to2; = —1 that corresponds 
to the original lA = lA-i. In this process, however, bounds on the rates of convergence are 
lost. This is the reason why we obtian in all of our results better bounds in the case of 
negative-definite functions h. 


3. Detailed proof of Theorem [3] 


Representation of correlation functions and truncation. 

We begin with a different representation of the correlation functions. For this we need 
a slightly generalized invariant ensemble: Define for a continuous Q of sufficient growth, 
continuous / of moderate growth and M G N the density on by 

1 


P. 


M 
N,Q,f 


(x) := 


n 


yM 

N,QJ l<i<j<7V 


=1 Q(^d+ETi/(^P 


X j X 7 


(29) 


We will usually have M = N+k for some k. If M = N, we will abbreviate Pn,qj ■= Pnq / 
^NQ ^NQf if / = 0- If / = 0 M = N, we have Pn,q- The A:-th correlation 


and Pj^ — t:>m 
function of Pn,v at points ti,... ,tk can be written as 

PN,viP^ • • •) 4 ) 

1 


I 

Jm 


N 


RN-k Zjsiy 


exp 




j=k+l 


i<j\ i,j>k 


X exp{ - N'^V{tj)+ 2 ^ log\ti-tj\} 

j=l i<j‘i 

xexp{2 log\ti — Xj\]dxk+i... dx]\f 

i<k, j>k 

= F{t) ^~^’^ Ej((_fcyexp{2 ^ log\ti-Xj\}, where 


Zn.v 


i<k, j>k 


F(f) :=exp{-iV^I/(tj)-F2 ^ log|fi-4|} (^0) 

j=l i<j\ i,j<k 

is the factor (f5n|) . which depends only on the fixed particles and y is the expectation 

of Pjj-k V label the random eigenvalues of the ensemble Pjj_^. y by Xk+i,... ,xn and 
abbreviate (ti,..., tk,Xk+i^... ,xn) by {t,x). Setting 


R{t,x) := R^_^ y{t,x) := 2 ^ log\ti - xj\ + log [P(t) 

i<k, j>k 


yN 

^N-k,V i 

Zn,v 


(31) 
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we arrive at the shorthand 


PN,vih, • • ■ , 4 ) = •)}• 


(32) 


Using 


and 


and recalling that we set U = 14 and U = Uu, we see that the fc-th 


correlation function of g at U ..., tfc can be written as 


Pn,q(^^i ■ ■ ■ i^k) — 




(33) 


Etv.u exp {W} 

The representation (I33p allows to control the effects of truncation on the correlation functions 
(see |GV14[ Lemma 6.3]). More precisely, let denote the expectation w.r.t. the ensem¬ 

ble P^y-i obtained by normalizing the ensemble P^y restricted to [—L,L]^. Furthermore, 
let Rl he the analogue of R in which the partition functions in (1311) are replaced by their 
truncated versions. Then |GV14l Lemma 6.3] (see also |Joh98j . [BdMPS^ l states that for 
each k there are L,C > 0 such that for all N and for all ti,..., 


Pn,q(^^^ • • •) 4 ) 


1 


^N-k,V-,L exp {(^ + -^l)(T •)} 


< e 


-CN 


(34) 


EAr,y;Lexp [U] 

Moreover, we learn from the last inequality in the proof of |GV141 Lemma 6.3] that L can be 
chosen such that in addition there exists a constant c > 0 with 

\Pn,q(^)\ < c (35) 

for all t G \ (—oo,L]^. It can be shown (see |GV14l Remark 4.5, Remark 4.8] and the 
proof of [GV141 Proposition 4.7]) that for such a choice of L and any A > 0 we can again find 
constants Gi(A),C' 2 (A) > 0 with 

0 < Cl (A) < EN,v-,Le^ < C2(A) (36) 

for all N, provided ag > a(A), where a(A) can be chosen as in (1241) and does therefore not 
depend on L. 

Linearization and proof of Theorem [3] b). 

We now give a more detailed description of the linearization procedure for negative-definite 
h. In this case we can indeed view —h as covariance function of a centered stationary Gaussian 
process on R such that (|26l) holds. Analyticity of the sample paths can be deduced from the 
spectral representation of / which we now explain in some detail. 

Let {Bt)t>o denote a standard 1-dimensional Brownian motion, in particular Bq = 0 
a.s.. Recall that by the law of the iterated logarithm, we know that \Bt\ is almost surely 
0{^/W\og^o^). Thus, for g : [0,oo) —)■ R of sufficient decay and smoothness, we can 
define the Wiener integral 


> roo 

g{s)dBs ■■= - Bsdg{s) = - 
Jo 


Bsg'{s)ds. (37) 

/o Jo Jo 

Recall that {Bt)t is a Gaussian process, that is, its finite dimensional distributions are Gauss¬ 
ian which is equivalent to the property that finite linear combinations of the family of random 
variables {Bt : t > 0} are Gaussian. By a limit argument we see that (|37l) is a Gaussian 
random variable and by elementary computations we find that the mean is 0 and the variance 
is /o°° 9{s)^ds. 
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Now, for the representation of /, let (-B/)t, denote two independent Brownian mo¬ 

tions. Then we can define 

f{t) := ^ J cos(ts)yJ—h{s)dBl + ^ J —h{s)dB‘^. ( 38 ) 

To verify that / defined in this way has the desired properties, it is enough to note that 
the right hand side of (j38p forms a Gaussian process on M (which can be easily checked 
using the characterization mentioned above) with mean 0 and covariance function —h. It is a 
somewhat surprising and maybe not so well-known fact |LS52] that the Fourier transform of a 
positive-definite real-analytic function decays exponentially at infinity (see also the discussion 
preceding (IbTT) !. By this exponential decay of h, representation (l38]l continues to hold for w 
from a strip {x + iy : x £ E., \y\ < <5} for some 6 > 0. For later use, we define the region 
D := (—L — 1, L -|- 1) X {—6/3, 6/3) C C with an appropriate choice of L made precise below. 
We thus see that / is analytic on D a.s.. It follows also from (f38]l that the extended process 
if{w)) is a complex-valued centered Gaussian process and it is straightforward to show 
that the covariance function is 

^{f{wi)f{w 2 )) = -h{wi - W^). ( 39 ) 

Before starting the proof we finally introduce an abbreviation for the Airy kernel determi¬ 
nants (see also (|7|)) 

Alkit) := det [A:Ai(ti, ij)]i<ij<fc • (40) 


Proof of TheoremlM b). Let k be fixed and choose L such that (IM)) and (f35]) are satisfied. 
We first restrict our attention to the case t £ [q, with 0 < e < min(2/3 — 2cj, 2/15). Note 
that this includes the region where the Airy kernel determinant describes the leading order 
behavior. Using (IMl) and the lower bound of (|36]) . it suffices to show 

{c*)-^N^/^E/!,_,y.^L exp {(A + RL){b + •)} - IEiv.y;L exp {U}AUt) (41) 

= O (N-^) 


uniform in t. Let / be defined as in ([381) and set f := f — f fdy. (cf. (f271l l. We now apply 
the linearization procedure and obtain from (1281) 


Observe that 


h,k /,\ 
Pn,Q]L\^) — 


^Ny-f/N-,LPN,V-f/N-,L^ 


^^Ny-f/N-y 


N 

Z^y-f/N-y = ^v,U;L IEa,U;L exp ( 


together with ([26|) leads to 


(42) 


EZ7v,y_//Ar;L = ZNyyENyyex.p{U). 
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The last three relations together with a truncated version of 

t 


yield 

h,k 


+ RL)ib + •)} = IEA,U;Lexp(ZY)p^Q.^(6 + 


N 


=E 


^N,V-,L exp { ^ / {Xj)}p%yj.L {b + ^^^2/3 ’ • • • > ^ + ^*^^ 2 / 3 ) 
i=i 


4 


where is the k-th correlation function of the unitary invariant ensemble PN,V-f /N-,L =• 

PNy,f \L defined on [—L, L]^. We thus get that (HTD is equal to 
N 

^Nyy exp { ^ /(xj)} (^(c*) ^N^^^P%yfy{b + ^^^2/3’ • • • > ^ + ^*^2/3) “ 


E 


i=i 


(43) 


By Proposition [16] (observe that c* = yv and ay = —by), by the determinantal relations 
and by the almost sure analyticity of /, we obtain for the term in the inner parenthesis 

{c*)-^N^/^p%yj.L[b + >•••>&+ ^*^ 2 / 3 ) = 0 {N-^) 


for almost all / with the O-term uniform for ||/||d < with k := 1/3 — a — ej^l. Since 
KKjyy-L exp { J2j=i — ^Nyy exp {W} is uniformly bounded in N (cf. (f36]l h we have 

derived 

N 


E 


i=i 


X 


^)->^N^/^P%yjy{b + 


h 


... ,b + 


tk 


r) - 


Hence the proof is finished for the case t G [q, if we can show 

N 


= 0{N-^). (44) 


E 


l{||/||o>V'^} X IEAry;iexp 


i=i 


X ((c*) ’'^^^^PNyf-,L{b + -^j^^,---,b+-jj^)-AIk{t)y -0{e ) (45) 

for some c > 0 uniformly in t G [q,N^]^. This bound will follow from an application of 
Holder’s inequality to separate the three /-dependent factors in (1451) . Here we will use 
norms for convenience. 

We first treat If follows readily from (l38|) that real and imaginary parts of / 

on D are (real-valued) centered Gaussian processes. By Borell’s inequality (see e.g. [ATOTl 
Theorem 2.1.1]) the supremum ||Al||oo of a real-valued continuous centered Gaussian pro¬ 
cess Xf over a compact K is sub-Gaussian, i.e. dominated by a Gaussian random vari¬ 
able with a certain expectation and variance a := sup^g^EX/. To apply this fact, we will 
without further notice always bound the supremum over the open set D by the supremum 
re sum of sub-Gaussian random variables is also sub-Gaussian, we see using 


over D. As t 


sup^eD 


fH 


< sup. 


'wGD 


Ref{w] 


+ sup^g£, 


Im/('u;) that sup^g^ f{w] 


has sub-Gaussian 
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tails. Since 


I fdfJ^ 


< sup^go 


f{w) we conclude that sup^g^j |/(t(;)| is sub-Gaussian. 


Note that |Re/| and |Im/| are both bounded by |/| and that K\f{w)\‘^ = —h{2ilmw) by 
(l3^ . We therefore find that sup^g£)|/(r(;)| is dominated by a Gaussian with a variance 
4sup^gj_5/3 5 / 3 ] —h{2iw). Thus there exists c > 0 such that 


P{\\f\\D > N'^} = 0{e and hence = C>(e (46) 

We now turn to the second factor Kn,v-l exp { ™ (IIS])- PropositionfT^ extends 

to the case of a truncated ensemble, yielding an error of exponential order which we omit in 
the following (cf. |GV14[ Remark 4.5]). Thus we have 


E 


N 2 

lEA,y;Lexp{ ^/(x/)} < Eexp ||[_l,l])}• 

i=i 


(47) 


Observe that the Lipschitz constant is evaluated over [—L, L] instead of M. 

In order to estimate this norm, we argue in addition to the sub-Gaussianity of ||/||[-l,l] 
that also the processes ||/dl[-L,L]) 11 /^^^II[-l,l] sub-Gaussian. To this end note that the 
process /' is a centered stationary Gaussian process with covariance function h". Borell’s in¬ 
equality yields again sub-Gaussianity of j/j^ and analogous arguments prove sub-Gaussianity 
of ll/^^^ll[-L, ^] as well. This shows the finiteness of the r.h.s. of ()17)l provided that aq and 
hence ay is sufficiently large. It is noteworthy that the condition on aq is determined by the 
variance h"{Qi) of f only. As f is stationary, this aq is independent of L and hence of k. In 
fact, this is the very reason that in Theorem [3] the convergence of all correlation functions 
can be derived for a given function h. 

Next, we will estimate p%yj.L{b+ ... ,b + -^y ^). From ([3]) we see that 

=■ A (where Kjy is now a shorthand for the kernel Kj^yjy associated 
to Pn,vj-,l evaluated at rescaled variables b + ’ 1 < J < ^) is positive semi-definite 

and can hence be written as A = for some Hermitian matrix B. Now using Hadamard’s 
inequality we obtain 


k k k 

det A = (det i?)^ < ^jj- 

j=ii=i j=i 


(48) 


In our case this reads 


Pkv,f;L(f>+ yJ^^2/3) - (^ - ^N(tj, tj) 

i=i 


- n + 

i=i 


G- 


y]\f2/3^’ 

Now, we use the well-known (see e.g. |Fre69] i representation 

,-NV+f ._ 


(49) 


„-NV+f 




















EDGE STATISTICS FOR REPULSIVE PARTICLE SYSTEMS 


18 


where the infimum is taken over all polynomials Pn-i of at most degree N — 1 with the 
property that P/v_i(t) = 1. From this it is obvious that 


Pn,vj-,l{^ + ^*^ 2 / 3 ) - PNy-A^ + ^*^ 2 / 3 )' 




By Proposition 1161 (|^. and by the uniform boundedness of the Airy kernel in the region of 
interest, we find that 




—-^T^) < and thus 


h 

c*Ny^ ’ ■ ■ 


■ ,b + 




(50) 


where C does not depend on /. 

This estimate together with the uniform boundedness of AI^ and the sub-Gaussianity 
proves the uniform boundedness of the third factor 


(c* 


-k]\^k/3 ^k 


Pn,vj-,l 


(6 + 


h 


(.j\r2/3 


,... ,b + 


tk 


^]\f2/3 


) - Alfc(t) 


in (I45p in the norm. This completes the case t G [q, 

Assume now that at least one component of t, say tj, is larger than N^. Then the 
asymptotics of the Airy kernel (see e.g. (I79p l imply that all entries of the j-th column 
of {Kp^i{ti,ti))i<i^i<k and consequently Alfc(t) are bounded by 0(exp(—A^^)) (cf. Remark 0] 
a)). We are left to prove a similar bound for the l.h.s. of ([9]). If one of the components tj is 
such that b+ ^,^ 2/3 > this follows from p35l) and for the remaining cases the estimate read¬ 
ily follows from Proposition [T7] a), (1491) and Proposition [16] for those components of t which 
lie in [q,N^]. Of course, Proposition [16] can only be applied for ||/|| < N^. One may use the 
same arguments as in the case t G [q,N^]^ to bound the contribution of {/ : ||/|| > N'^}. □ 


Extension to general h. 

If —h is not positive-definite, we may write it as a difference of positive-definite functions. 
Denoting by g± nonnegative and negative part of a function, we first write h = (/i)+ — (h)_. 
Setting as the inverse Fourier transform of h±, we get a decomposition h = — h~ oi 

h into positive-definite, real-analytic functions. It is this step where the assumption of expo¬ 
nential decay of h is needed. Exponential decay of h is equivalent to being real-analytic. 
Sufficiency is easily seen as the exponential decay allows for an entension of the Fourier repre¬ 
sentation of h to a strip D from which analyticity of can be deduced. For the necessity we 
remark that any real-analytic positive-definite function has an exponentially decaying Fourier 
transform [LS52[ Theorem 2] and thus with also h must have this property. 

Define for a complex parameter 2 ; G C 

N 

=|( “ Kk]) (hi) 

+ 2 ^ ^ (xi - Xj) - [h^ {xi) + {xj) - ). 


(52) 
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We have lA-i = Li. As seen in (1411) . we have to show 

exp{(pi, + RL){b + •)} - Ea,U;L exp {U,}Ahit) = o(l) (53) 

for z = —1 uniformly in t E [q,oo)^. Let us recall two basic facts from analysis. As our 
linearization procedure only allows for nonnegative real z, we will prove (j53|l for positive real 
z and use complex analysis to deduce (fM]l also for z = —1. Recall Vitali’s Theorem for 
instance from |Tit39( 5.21]: Let {fn)n be a sequence of analytic functions on a domain U C C 
with \ fn{z)\ < M for all n and all z £ U. Assume that lim„^oo/n(^) exists for a set of z 
having a limit point in U. Then lim^^oo fniz) exists for all z E C/ and the limit is an analytic 
function in z. 

To capture the uniformity of the convergence in (I53p in a way which preserves analyticity 
in we use the following obvious characterization. A sequence of complex-valued functions 
(/n)n converges uniformly on the sequence of sets An C towards a function / if 

and only if for all sequences {nm)m C N with lim^^oo nm = oo and all sequences {tm)m with 
tm £ -^rim have limm—)-oo fnmi.^m} f itm) — 0. 

For our application, we take A := A^v := [q,oo)^. Let {Nm)m C N be a sequence going to 
infinity and {t(rn))m be a sequence with t(^nn) ^ Define Wm '■ C — C by 

W„(Z) ■McT'‘lVi/^KZ-Z.V:L + Rl)(I> + •)} 

C lym 

-'^Nm,V;Lexp{Uz}AIkit^jn))- ( 54 ) 

With these two observations, we are ready to complete the proof of Theorem [3l 


Proof of Theorem\^ a). We start by showing convergence of Wm{z) uniformly for z E [0,1]. 
We will mostly omit the index m in the following. Define —h^ ■= zh~^ + h~ such that /i_i = h, 
denote by fz the corresponding Gaussian process, and set fz ■= fz — J fzdjj,. Now, using the 
same arguments given in the proof of Theorem [3] b) we find that 


N 


E En,v-l exp {y^Jz 
i=i 

= 0{N-^). 


1= ^j—k j\^k/3 ^k 


PNyj; 




tjY2/3 


... , 6 - 1 - 


ik 




) -Al, 


The convergence is uniform in 2 E [0,1]. This is a consequence of the boundedness of 
suPuje[-(5/34/3] ~hz{‘^iw) for 2 ; E [0,1] (cf. the arguments above (j46]l i. Note that the choice 
of aq in the hypothesis of Theorem [3] depends on this bound. This settles the convergence 
of Wm on [0,1]. 

We complete the proof by demonstrating uniform boundedness of Wm on the domain 
G := {z £ C : Rez < 1}. We first observe that applying ([2X1) and (l2^ to /i+ instead of h yields 
the non-negativity of ~ ^j) ~ ~ ■ Thus ReUz{x) < Ui{x) 

for all 2 ; E G and all x and consequently we have 

supE7v,y;L |exp(6/^)] < supEat^;!, exp(6/i) < 00 , (55) 

N N 

where the last inequality follows by linearization as above. Hence the second term in the 
definition of Wm is uniformly bounded on G. Furthermore, the uniform boundedness of the 
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first term in (I54p follows from < lAi for all z ^ G, from the convergence and hence 

boundedness of VFm(l) and from We are now in a position to apply Vitali’s theorem, 

providing the desired convergence of VFm(—1) to 0. □ 

Finally, we would like to point out a subtlety that can be explained e.g. by looking at 
relation (133I1 . The z-dependence used above is introduced by twice replacing U by Uz on the 
right hand side of (|33p . Observe that this is in general not the same as replacing h by hz 
on the left hand side. In fact, we never consider the ensembles with z ^ —1, i.e. for 

hz 7 ^ h. The function V and the associated equilibrium measure fj, do not depend on z either. 

4. Proofs of Results on the Largest Particle 
Proof of Theorem O We have to compute the limit of the gap probability 

PN,Q{xj ^{b+ J = 

This proof uses the same techniques and route as the proof of Theorem [3l However, in that 
proof the truncation threshold L depends on the order of the correlation function k, leading 
to an O-term depending on A: in a non-obvious way, e.g. in the expectation of quantities like 
WfWn or \f\jr- It is therefore convenient to truncate the event before expanding in terms 
of correlation functions. To this end, we recall from the last inequality in the proof of the 
truncation lemma [GV141 Lemma 6.3] that 

PN^Qi^) ^ expjCA^ - ciN[V{t) - C 2 log(l + t^)]} (56) 

for some positive C, ci,C 2 . Hence, if L is chosen large enough, we have for some c > 0 

POO 

PN,Q{xj i [-L,L] for some j) < 2N = 0{e~‘'^). 

From this we also conclude that the replacement of the normalizing constant q by q.^ 
is negligible, where Z^q.^ denotes the normalizing constant for the ensemble P^q.^ that 
arises from Pjj-g by restricting it onto [—L,L\^. We thus have 

PN,Q{xj i (6+ -^^,oo) for all j) = P^,Q.L{xj i (6+ for all j) +0(e-"^). 

(57) 

To evaluate the latter probability, we proceed as in the proofs of Theorem [3] a) and b). We 
will write explicitly only the case of negative-definite h, the general case follows then with the 
additional arguments given in the proof of Theorem [3] a). Note that the seemingly stronger 
uniformity for s G M in a) is a simple consequence of the continuity of F 2 . In order to show 

Pn,q-A^j i (6+ -^,L), j = 1, ...,iV) - F2{s) = 0{N-^) 

it suffices to prove the bound 
N 

E [EMy;Lew{Y.f{xj)}{PN,vj;LA i (6+ -^,L),j = 1,..., N) - ^ 2 ( 5 ))] = 0{N-<^) 
j=i 

(58) 
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(cf. In order to see this, integrate over [—L, and use ()12]) together with 

the lower bound of (1361) . As in the proof of Theorem [3] b) we choose e and k depending on cr. 
Again we can neglect the contribution of those / to the expectation in (j58p with ||/||_d > N^. 
We are therefore left to prove 


PN,v,f;L{xj i (6 + -^,L),j = 1,...,N)-F2{s) = 0{N-<^) (59) 

uniformly for all / G Xo with ||/|| < N'^ and all s G [q, {L — h)c*N'^P]. To this end, we will 
represent both probabilities as Fredholm determinants. It is well-known (see e.g. |TW98] for 
a nice derivation) that 

PN,v,f-,L{xj ^ {b + = det{I - }CN,v,f-,L) =■ 

where JCN,v,f;L denotes the integral operator on L^((s, (L — b)c*N^P)) with kernel 

Kn,vj-,l{PP 2) := -j;^^PN,v,f-,Lib + ^^^,b + 


Observe that K agrees with the definition in Proposition [16] since b = by and c* = 71 /. For F 2 
we have F 2 {s) = A(/CAi) with /Cai (cf. (fTU]) i being the integral operator w.r.t. the Airy kernel 
on L?‘{{s, 00 )). For comparison, it will turn out convenient to consider Fredholm determinants 
of integral operators that are defined on the same L^-space. To this end, define KLai^l as the 
integral operator on {L — b)c*N'^P)) w.r.t. the Airy kernel. Let us show that replacing 

P 2 {s) by A(/CAi;L) only results in a negligible error. Let C, denote the determinantal point 
process on M determined by the Airy kernel (see e.g. [AGZlOl Proposition 4.2.30]). Then both 
F 2 {s) and X{K.a\-l) can be represented as gap probabilities w.r.t. C- Let P be the probability 
measure underlying C,. From the asymptotics (|llh of F 2 we conclude 


|F2(s)-A(/CAi;L)| 


P(c n (s, 00 ) = 0) - P(C n (s, (L - b)c*NP^) = 0) 


< P(Cn [(L - b)c*NP^,oo) / 0) = 1 - F 2 {{L - b)c*NP^) = ©(e"^^) 


for some c > 0 . Thus it suffices to estimate the difference A(/Ca,i/,/;l) — A(/CAi;L)- Here we 
use again the classic series representation of Fredholm determinants [AGZlOl Def. 3.4.3] and 
the inequality [AGZlOl Lemma 3.4.5] on the difference between the Fredholm determinants 
of two integral operators S, Q with bounded kernels S,Q on an interval (c, d). Let u be a 
hnite measure on (c, d) with total mass IJujji. Let A(5, u) := det(/ —5) denote the Fredholm 
determinant of S on L?‘{{c, d), v). Then 


lA(5,u)-A(Q,u)l < 



^i+fc/2||^||fcmax(115’lloo, llQlloo)^ 

k\ j 


ll^-Qlloc, 


(60) 


where jj • jjoo denotes the sup-norm on (c, d)^. A natural choice would be to use the Lebesgue 
measure for v on the intervals (s, {L — b)c*N‘^P). However, Uujji would then be of order 
and consequently the series in (I60|) could not be bounded uniformly in N. For our application 
one can circumvent this issue by transfering the fast decay of the kernels onto the measure 
ly. In order to do so we work with Lp‘{{s, {L — b)c*N‘^P), u), dv{t) := e~‘^^dt and set 

SN{tl,t2) ■= , Q{tl,t2) '■= KAi{ti,t2)e^^~^^^ . 
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Using the above mentioned representation of the Fredholm determinant one immediately 
obtains /\{Kn,vj-l) = ^{^Nyj-,L,dt) = A{SNy) and A(/CAi;L)= A{)CAi-L,dt) = A{Q,u). 
Observe that ||p||i is now bounded uniformly in N. The same holds for ||Q||oo which follows 
from the asymptotic behavior of the Airy kernel as presented e.g. in |KSSV14l (4.23)]). The 
uniform boundedness of USArjloo can be derived as follows. Note first from the determinantal 
formula (jl]) for k = 2, from the positivity of the 2-point correlation function and from the 
symmetry of the kernel that 

*S'Ar(UT 2 )^ < SN{tl,ti)SN{t2,t2). (61) 

The boundedness of Sw(t, t), t E (s, {L—b)c*N‘^/^) on the diagonal can in turn be derived from 
the first two statements of Proposition [16] and from Proposition [T7| a), where the asymptotic 
behavior of the Airy kernel and a lower bound on rjv{x) of the form c'(x — 1)^/^ are being used 
in addition. We have now shown that the series on the right hand side of (|60p is uniformly 
bounded in N if we choose S^i for S. Establishing (1591) uniformly for s E [q, (L — b)c*N‘^^^] 
thus reduces to proving IjS’^r — Q||cxd = To this end recall e = 2(| — a — k) > 0. 

We distinguish two cases. For ti, t 2 E [o', Proposition [16] yields the desired bound. In the 
case that one of the variables is larger than N^, the crude estimate I^at — Qj < \Sn\ + |Q| 
together with ()61l) and Proposition [T7] a) provides a bound of the form that 

suffices easily. □ 


We continue with the proof of our result for the regimes of moderate and large deviations 
for the upper tail of the distribution of the largest particle, which combines the analysis 
devised in [Schl5[IEKS15j and the asymptotic results of |KSSV14] with the procedure of the 
previous proof. 


Proof of Theorem 0 We start by observing that it suffices to prove the statements of Theorem 
[ 6 ]for N > Nq and in the cases of statements a), c) for t E {b + CN~‘^^^,T) for arbitrarily large 
but fixed constants Nq and C. Let us first turn our attention to statement c). We begin with 
the truncation of the ensemble to [—L,L] as in (I57p . It follows from (|56[) that the constant 
c > 0 in the error 0{e~'^^) of that truncation can be chosen arbitrarily large by increasing 
L. Thus, we may choose L > T + 1 so large (depending on the choice of T in the statement 
of the theorem) that we have for all t E (6 -|- ,T)-. 


^N,Q (^max > t) ^ IEAr,V;L [lA{x)e^^^'>] Ii i CO (' -cN\] 
T'Nyit) ENy-L^^^^J^Nyit) 


(62) 


for some new positive constant c and with A := {x € : rcmax > 

Now we apply the stochastic linearization to rewrite the numerator and obtain (use in 
particular the relation for the partition functions provided by (|42p l 


Eat.i 


1 Jrr\MhA 


In view of ([ 6 ^ , (f63ll we first derive upper and lower estimates on PNyjyixmia > t)/FN, y{t)- 
Since PNyj-y is determinantal we have 


PNyjyixmsix > t) — 1 — A{}C]s[y,f-L), 


( 64 ) 
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where A(JC]sfyj-L) denotes the Fredholm determinant related to the integral operator with 
kernel and acting on L), dx). Next we apply Proposition [T71 bl. The formula 

there holds for all t G {b + , T) (since we have ensured above that L > r + 1) and for 

all / G Xd satisfying ||/|| = ||/||_d < coN(t — b), where cq, Cq denote some positive constants. 
In this regime a first simple consequence is 


KN,v,f-,L{y,y)dy = O 


1 


N\t-bf‘^) ’ 

where one should use from the proof of Proposition [T7] b) that 

^-Nr,v(t/bv)+0{Vt^\\f\\) ^ < 1 . 


The definition of the Fredholm determinant as in [AGZ10( Def. 3.4.3], Hadamard’s inequality 
(see ((48]) ) and Proposition [T3b) then yield the finer estimates 


1 - A{ICNyf-L) = KNyfy{y,y)dy^ 


1 + 0 


1 + 0 


N\t-bf^^J 

+ 0 


N\t-b\ 


^N\t-bf^^ 


Recalling (f 6 ^ we have derived for all f G (& + CqN the upper bound 

-P/V,U,/;L(^max ^ ^ , Ol i-in 

-- <1 + 7(11/11), With 


(65) 


( 66 ) 


7(r) 


c 

Cr^yt — b 

gCr+t—b 

1 

, -+v.v(b 


for 

for 

for 

for 


iV |t- —0| 


< r < 


N\t-b\‘‘ — ' ^ Vf-b 

^ < r < coN{t - 
coN{t — h)<r 


(67) 


for some suitable constant C. In the case r > CQN{t — b) we have used that the numerator of 
the left hand side of (f 66 |) is a probability and thus bounded by 1 . 

We now formulate a basic integration lemma that is convenient for estimating the right 
hand side of (j63p . To this end we introduce a function F that is related to the distribution 
function for ||/|| induced by the Gaussian measure on / and modified by the positive density 


F{a) := E _ 


( 68 ) 


Lemma 15. Let 0 < a < (3 < oo and assume that 7 : [a, /3) —>■ [0, 00 ) is a continuously 
differentiable function. Then 

+ 7 ( 11 / 11 ))) < F{a) - F{/3) + 'y{a)F{a) + J F{r)ff{r)dr. 

(69) 
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To see this, observe first that the left hand side of (f69l) is given by the Riemann-Stieltjes 
integral 


rP r0 

/ 1 + 7 (r)(iF(r) =-(1 + 7 )F 1^ + / F{r)'y'{r)dr 

J OL J a 


after integration by parts. Using —'y(l3)F{l3) < 0 completes the derivation of (l 6 ^ . 

We apply Lemma [TSl on four intervals [a,/3) = [aj-i,aj) for j = 1 , 2 , 3, 4 according to 
the cases in (|67|l . i.e. oq := 0, ai := «2 := 7 ^) ■= coA^(t — b), and 04 := 00 . 

Combining ([62|) . (I63]l . (I66l) . (f67)l and using = F(0) ( SG 6 (j26p ) WG obtOjin 


Ejv,U;L [lA(x)e^(")] 


< 






We are left to show that the last six summands on the right hand side are of the form 
O (^1/{N |t — or of the form 0{^/t — b). Since 7 ( 0 ) and 7 ( 01 ) are of the first form and 

trivially F{ai)/F{Qi) < 1 the first two of these summands are good. Since F{r)dr < F{0) 
it is clear that the hrst integral in the sum is of the second form. For the remaining three 
terms we use that there exist positive constants d, K such that T(r)/T(0) < Ke~'^^ for 
all r > 0 and for all N. To see this, use the sub-Gaussianity of ||/||d analogously to the 
derivation of (|45]) together with the lower bound on F{0) = K]\fyy^ provided by (f36]l . 
Then it follows 


7(a2)T(a2)/T(0) = O = 0{Vt^). 


Moreover, setting x := N \t — bf^'^ > Cq > 1, y ■= x^!"^ and ci := cZcq, we have 


for N sufficiently large (e.g. > 4/ci would do) and is therefore of the first form. Finally, 

we consider the last integral in the sum. It can be estimated by 

O (^Vt^ J = O (yt^/{da 2 )^ = O • 

This completes the analysis of the upper bound. The lower bound is proved in exactly the 
same way. One only needs to replace 1 + 7 by 1 — 7 . In summary we obtain 


Eiv,U;L [lA(x)e^W] 


l + O 


\N\t-b\ 




(70) 


for all t G (6 + ,T). In view of the remark at the very beginning of the proof and 

using ([62]) we have established statement c). 

We turn our attention to the case of general functions h that are not necessarily negative- 
definite. Observe first that ([70t) still holds if U is replaced by Uz for any positive z where we 
do not insist on uniformity of the O-terms in z. To see this one may use exactly the same 
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proof of as above with F being replaced Fz which is defined via (l 68 ]l with the Gaussian 
process fz instead of /. One may again show that the tails of Fz are sub-Gaussian and that 
infArF;j(0) > 0. These were the crucial properties of F used in the derivation of (fTOj) . Using 
in addition (l55]) we conclude 


IEjv,y;L [lA(x)e^-(")] 


—>-0 for ^ oo 


(71) 


for 0 < z < 1 and for those values of t that are considered in statement b) of the theorem. 
For later use we note that by the same reasoning, the left hand side of (I7ip remains bounded 
for 0 < z < 1 and all t G (6 + CoN~^^^, T). In order to apply Vitali’s theorem and in view of 
(|55]) it suffices to show boundedness of 


^N,V;L [lA(x)e^-(-)] 

e.g. on G = {z G C : Rez < 1}. Using the pointwise estimate Rel/(z{x) < Ui{x) for all 2 ; G G 
this follows from (EH) for z = 1 and from (I55|) . Vital! now yields convergence EH) also for 
z = —1, i.e. for Uz being replaced by U. The lower bound of (IH6]l for A = 1 together with 
(162]) prove statement b). 

Observe that we have already shown the upper bound of a). In fact, by (l6^ and by the 
lower bound of (j36n for A = 1 we are only required to prove an upper bound on (1721) for 
z = —1. This can be achieved as above where one needs to use the boundedness of m on 
all of (6 + GoV-2/3,r). 

We establish the lower bound by contradiction. Assume that there is no such bound. Then, 
again by (f62]) . for any fixed (large) G > 1 there is a sequence {t]s[)N in {b + ,T) such 

that 


,. '^N,V-,L |1Ajv(®)' 

hm - :r 7 r-^ - ^—7 

Af^oo Em vLf^^^^J^Nv{tN) 


'i 


(73) 


where Am := {x G : Xmax > By monotonicity of Uz{x) for real values of z we have 
for 2 G (— 00 , —1] 


^N,v-,L [lA^(a;)e^W] ^ Em,v-l [lA^(x)e^-(")] 

J^Nyitn) ~ J^NyitN) 

and thus, by the boundedness of EMy-L^, both sides converge to 0. By Vitali’s theorem, 
this convergence extends to z G (0,1) as well. However, we will prove that for z = 1/2, (1721) 
is bounded away from 0 uniformly in N and t G (6 + ,T) for sufficiently large values 

of G, thus providing the desired contradiction. 

Recall the definition of -^ 1/2 below (I70p and choose /3 such that Fi^ 2 (/d) — '^i/ 2 (b)/ 2 . In 
view of (fHH) and (l65]) there exists G > 0 such that we have for all t G (5 + T): 


d{T) : = 


inf 

ll/ll</3 


PNy,f\L (^^max 


> t)/FMy{t) > 0 . 


Thus the right hand side of (f 6 ^ with U replaced by Z 7 i /2 and devided by FMy{t) is bounded 
below by d{T){l — Ti/ 2 (/ 3 )) > d{T)Fij 2 {R)/2. Since in addition ^ 1 / 2 ( 0 ) > F{0) = EMyye^ 
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is bounded away from 0 we have arrived at the advertised lower bound on (j72|) for 2 : = 1/2, 
completing the proof of Theorem [6l □ 


Proof of Corollary \11[ The corollary follows from an application of the Borel-Cantelli lemma. 
We thus have to show for any fixed e > 0 

00 

^A,Q(kmax - 6| > e) < OO. 

N=1 


The necessary estimate for the probability of Xmax > 6 + e is provided by Theorem [6] a). It 
remains to derive a bound for the lower tail of the distribution of Xmax- This can be done as 
follows. Given e > 0, let g he a smooth, nonnegative Lipschitz function with gdg = 1 and 
with support [b — e + 6,b — 6], where 0 < 5 < e/2. Then 




N. 


N N 

= 0 ) ^ Pn.q (; 

i=i i=i 


,Q(xmax <b-e) < Pjf^Q ( g{xj) = 0) < Pff^Q ( Y g{xj) -N / gdg>N). 


With (|2^ . (|23|) . and (l2^ an application of Holder’s inequality gives 



for some C > 0 and some A > 1. An application of Chebyshev’s inequality to the statement 
of Proposition [13] with a good choice of e shows that this last probability is of order 
for some c > 0, provided that N is large enough. The corollary is proved. □ 


Appendix A. Asymptotics for determinantal ensembles 


In this appendix we extract those results from |KSSV14] on the asymptotics of the de¬ 
terminantal ensembles Pn,vJ\L '■= PN,v-f/ n-,l^ defined on [—L, that are relevant in our 
context. 

We begin by recalling the notion of the equilibrium measure with respect to a general 
convex external field V. By this we mean the unique Borel probability measure /ip which 
minimizes the energy functional 


/i I—)• 


log|t —s| ^ dg{t)dfj,{s) + 


V {f)dg{t). 


A general reference on equilibrium measures with external fields is |ST97j . It is known that 
under mild assumptions such a unique minimizer exists. In our case V will be convex and 
the density of /xp can be described as follows (see e.g. |KSSV14] Sections 1 and 2]): There 
exist real numbers ap < 6p that are uniquely determined by the two equations for a and b, 



V'jt) 

■s/{b- t){t - a) 


dt = 0 , 



y/{b- t){t - a) 


dt = 271. 


We denote by Ap the linear rescaling that maps the interval [—1,1] onto [ap, by], 




et-c/ 




V 


t + 


V 


+ b-T/ 


V 


(74) 
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and introduce a function Gy on M by 


Gy{t) 



{V o XyY'it + u{s - t)) 

Vi — 


duds , 


(75) 


that inherits real analyticity and postivity from V". The density of the equilibrium measure 
is then given by (see [KSSV14( (2.1), (1.11-1.14)]) 

d'lJ'yit) = 77 ^ \l ~ ~ GjGy{\~ (t))l[(j^ (t)(it. 

[Oy tty) TT V 

The linearization technique of [GV14| requires to consider ensembles PN,v-f/ n-,l that arise 
from restricting the Pjv,y-//v to [—L,L]^ and renormalizing them as probability measures. 
Note that the choice of L throughout the paper ensures that the support [ay_f/]y,by_f/jy] 
of is contained in the interior of [—L,L]. Moreover, the functions / are always 

defined on some fixed domain D in the complex plane that contains [—L,L]. They belong 
to the real Banach space Xd := {/ : D — C : / analytic and bounded, f{D n M) C M}, 
equipped with the norm ||/|| := ||/||d := l/(^)l- |KSSV14l Lemma 2.4], the maps 

V ^ Uy, V ^ by that relate the external fields to the endpoints of the support of their 
equilibrium measures are G^ with bounded derivatives on a neighborhood of 17 in Thus, 
for sufficiently large values of iV, 


avj - av, bv,f - by = 0 , (76) 

where the subscript yj is short for y_fi]y. We can now formulate a first result on kernels 
Kn,vj-,l associated to determinantal ensembles of the form PN,v,f-,L that has essentially been 
proved in [KSSV14j . 


Proposition 16 (cf. [KSSV141 Theorem 1.8]). Let V : [—L,L] —)• M be real-analytic with 
inf^g[_£ £] 17"(f) > 0. Assume that the support [ay, by] of the equilibrium measure fiy is 
contained in {—L, L) and recall the notations introduced in this appendix. Moreover, define 


where s, t G M. Then we have for fixed 0 < k < 1/3, 0 < e < min(2/15,2/3 — 2 k ), (7 G M, 


PN,VJ-,L{s,t) 


'KM{s,t)+0{N--^P) , ifq<s,t<2, 

< KAi{s,t){l + 0{N^+P^-P^)) , ifl<s,t<Nfi 
KAi{s,t)+Af{t)0 , ifq<s<l-,2<t<N^ 


with the O terms being uniform in N, in f G X^) with ||/||d < o.nd in s,t within the 
respective regions. 


Proof. The result in |KSSV14l Theorem 1.8] provides leading order information on 


[by, + Kj + 
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together with error bounds that have the uniformity needed to prove Proposition [T6l We 
note in passing that (f76]l ensures that, uniformly for f £ Xd with ||/||_d < [ayj,6yj] C 
(—L,L) for N large enough. 

The reader should be aware of some differences in notation, e.g. the quantities Q, V, 7 y, 
Ai in |KSSV14] correspond to V, V — f /N, 'JvjiWj ~ ®u,/)/2, A'ai in the present paper. 
Our first task is to replace byj by by and yyj by yy. Using 

[(U - f/N) o \yjf - [U o Ay]" = [{>!yj?V'' o Ayj - (A'y)V" o Ay] - {X'yjf{f''/N) o Xyj, 

relations (|74p - (|76p . the thrice differentiability of V on [—L,L], the analyticity of / in D, and 
the positivity of Gy, one obtains 

Gyj{t) = Gy{t){l + 0{\\f\\/N)) (77) 

uniformly for f £ Xd with \\f\\D < 0 < ac < 1 and also uniformly in t E M with both 

Ay(t), Xyj{t) £ [—L,L]. This yields in particular yyj = 7 y(l + 0{\\f\\/N)). In order to be 
able to apply the results of |KSSV14] we relate to any q <t < a number t = t{N, V, f) by 

+V2757;7 

and in the same way we relate S' to s. For all relevant values of /, t, and s we obtain the 
uniform estimates 

t = t + tO{\\f\\/N) + N‘^/^-fyj{by - byj) = t + ^ s = s + (78) 

Next we discuss the different regions in the {s,t) plane one by one. 

(a) q < s,t <2: Application of Theorem 1.8. in |KSSA714j yields 

KNy,f-,L{s,t) = ^ (i^Ai(s,?) +0(fV-2/3)) 

and the claim follows from (|78p and from yyj/yy = 1 + Observe that we need a 

slightly modified version of the result in |KSSV14] , because one only has that 's, t lie in some 
small neighborhood of [q,2]. However, the proof in |KSSV14] shows that the statement also 
holds true in such an slightly enlarged neighborhood. This remark applies to the next two 
cases as well. 

{h)l<s,t < N^: Theorem 1.8 of |KSSV14] together with e < 2/15 show 

KNy,f-,L{s,t) = KAi{s,t) (l + _ 

The asymptotic formula for the Airy kernel (see e.g. [KSSVl^ (4.23)]) implies the existence 
of some G > 1 with 

^ < XAiis, t)(st)^/^(x/s + Vt) exp[|(s^/^ + t^/^)] < G (79) 

for all s, t > 1. Using e.g. |KSSV14l Proposition 4.2] together with the asymptotics of the 
Airy function and its derivative (see |AS64[ 10.4.59, 10.4.61]) one obtains a similar formula 
for the partial derivatives of the Airy kernel leading up to 

\VKAi{s,t)\=KAi{s,t)0 (y~s + Vt^ =KA,{s,t)0 
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Since e < 2/3 — 2k relation (17^ implies for all s*, t* with distance 0{N^ from s, resp. t 

KM{s*,t*) = KM{s,t)0{l), 


yielding the claim. 

(c) q<s<l]2<t< N^: This time application of [KSSV14] gives for e < 2/15 

KNyj-,Lis,t) = [KmCsS) + (l + 0(iV-i)) . 

The standard techniques used to derive the asymptotics of the Airy kernel in [KSSV14] (see 
in particular Proposition 4.2 there) yield the following estimates for the values of s and t 
considered in the present case: 

KM{s,t) = Ai'(t)0(t-i), |Vi^Ai(s,t)| = Ai'(t)0(t-i/2), Ai'(f) = Ai'(t)0(l), 


where we have used e < 2/3 — 2 a to derive the last relation. This completes the proof. □ 


In the regime of moderate and large deviations Proposition [16] does not provide the leading 
order description, which is the content of the following proposition. As it turns out, it suffices 
to consider the Christoffel-Darboux kernel on the diagonal. In order to simplify the resulting 
formulas we now make use of the assumed evenness of Q and h that implies the evenness of 
V (cf. Remark [2] b)) and hence ay = —by- 


Proposition 17 (cf. [KSSYld] Theorem 1.5]). Let V satisfy the assumptions of Proposition 
[T^and assume in addition ay = —by. In addition to the notations introdueed in this appendix, 
reeall the definition of the funetion r]y in (USD. There exist constants cq > 0 and Cq > 1 only 
depending on V such that for 

a) by + CqN-^/^ <t<L: 


KNyjyit.t) = 


lj^f.-Nvv{t/bv)+o(Vt:=W\\f\\) 


h) by + CQN-^l^<t<L-l: 


1 + 0 


N\t-by\ 


+ 0 


N\t-byf^. 


J^I<:Nyj.L{y,y)dy = 


5^g-iVw(+fev)+C’(vTTP||/||) 
-b‘^)r]'y{t/by) 


1 + 0 


N\t-by\ 


+ 0 


N\t-by\^/^^ 


All O terms appearing in statements a) or b) are uniform in N, in t G {by + CqN 
resp. int£ {by + L — 1), and in f ^ Xr, with ||/|| = ||/||_d < coN{t — by). 


Proof. Prom the arguments given in the beginning of the proof of Proposition[T6|it follows that 
we can apply the result [KSSV141 Theorem 1.5 (ii)] to the ensembles Pny-f /n-,l provided 
that ll/ll/Y is sufficiently small which we can always achieve by choosing the constant cq in 
the statement of Proposition [T7| small enough. Hence there exists a positive number C such 
that for all X > 1 + with Xyj{x) < L 


KNyj-,L{^vj{x),\yj{x)) = 


27r{byj — ayj) 


„-Nr)vj{x)_ 


x-^ — I 


1 + 0 


1 


Y lx — 1 


3/2 
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holds with the required uniformity. Setting x = Ayj(t) we immediately arrive at 


K 


N,V,f-,L 




271 


- W,f){t - avj) 


1 + 0 


N\t-hvj\ 


3/2 


(80) 


Using again the smallness of ||/||/[A^(t — by)]-, we learn from (|76p that for real exponents a 
one has 


{t - hyjY = {t- bvY 


1 + Oo 


N\t-hv\ 


(81) 


and consequently 

4:{t-bvj){t-avj) 


1 + 0 


N^t — byj 


■372 


= fey 


1 + 0 


Af|i-bv 


+ 0 




(82) 

We now turn to the exponential term. Applying (1771) to the definition of r]v in (I15p . we see 
r]vj = r]v{l + 0{\\f\\/N)) uniformly on the domain of interest. Thus 

vv,f - w = w 

Moreover, it is straightforward to see that Ay^(t) = t/by, A(/^(f) — Ayj(t) = 0(||/||/iV), 
r]y{x) = 0{\x — 1|^^^), and 7?(/(x) = 0(|x — 1|^^^). Using in addition (fSTT) we obtain 

w,/ ( 1 ^ “ 

W ~ W (V^u) =C’ 

^-Nrivj{\Yjit)) _ ^-Nrjv{t/bv)+0(Vt-bv\\f\\) (-§ 3 ^ 


that leads to 


Combining (fHni) with (l8^ and (l83p completes the proof of statement a). 
The essential relation for establishing the second claim is 

p-Nr](t) 


I 


L Q-Nr)(y) 

't {y-b){y-a) 


dy = 


N{t — b){t — a)ri'{t) 


1 + 0 


1 


N\t-bf^ 


(84) 


with an O term that is uniform in N, t, and / in their respective domains and where a = ayj, 
b = byj and rj = rjyj o Ay^ We first sketch the derivation of (fMp following [SchlSi Lemma 
4.8] (cf. [EKS15| L Since rj is strictly increasing the substitution u{y) := rj{y)—rj{t) is invertible 
and we obtain 


I 


L Q-Nr]{y) 

{y - b){y - a) 


dy = / e-^^k{u)du, k{u) := 

Jo 


1 


(y(u) - b)(y(u} - a)r/'(y(u)) 

Observe that k(u) = k(0) + uk'(u) for some u G (0,tt). In addition, straightforward estimates 
yield 


k'(ti) 

k(u) 


= O 


1 


Jy(u)-bY/\ 


and by monotonicity of k and y: k'{u) = /c(0)O 


1 


}t-b\ 


3/2 

































We arrive at 
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o-Nri{y) 


It {y-b){y-a) 


dy = e-^^WA:(0) 




\t-b\ 


3/2 


with d = r]{L) — r]{L — 1) > 0. Evaluating A:(0) provides (l8ll) . 

Note that (I84p implies in particular that K]\fyj.L{y,y)dy has exactly the same repre¬ 
sentation as K]sfyj-L{t,t) in ([80]) except for additional factors N and (f/u,/ o the 

denominator. Thus the arguments in the proof of statement a) below (|8UI) together with 


ivvj ° - V,/)(t - av,f)Gv,f{Xyj 


it)) 


= -j^Vvit/bv) 


1 + 0 


N\t-bv\ 


prove statement b). 


(85) 
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